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ANNOTATION 


The unsupported phase of motion of a multileg bounding vehicle Is 
examined In the work. An algorithm for stabilization of the angular 
motion of the vehicle housing by change of the motion of the legs 
during flight Is constructed. The results of mathematical modeling of 
the stabilization process by computer are presented. 
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ALGORITHM FOR THIS STABILISATION OF MOTION A BOUNDING 
VEHICLK IN THIS FLIGHT PHAGIS 

V . V « Lapshin 

Institute of Applied Mathematics :1m. Keldysh, 
USSR Academy of Sciences, Moscow 


Introduction 


The problem of stabilisation of the motion of a bounding vehicle /G* 
in the flight phase is investigated in the work. An algorithm for 
stabilisation Ls a basic element of the motion control system of a 
bounding vehicle. 

A vehicle ls considered, which consists of a housing and four or 
six two section logs, each of which has three degrees of freedom. The 
total weight of the legs is a significant portion of the weight of the 
housing . 

The motion of the vehicle consists of alternation of two phases: 
support, during which all the legs stand on the supporting surface and 
there is quasistatic stability, and the unsupported or flight phase. 

In [1, 2] , a mathematical model of the three dimensional motion of the 
vehicle is constructed for the supported and unsupported phases of 
motion. In the flight phase, the first integrals of the equations of 
motion are obtained (the motion of the center of mass of the vehicle 
along a ballistic trajectory and the law of conservation of the angular 
momentum of the vehicle about the center of mass). 

In the support phase of motion of a bounding vehicle, high speeds 
and accelerations develop, for which great forces and power are re- 
quired. Therefore, it is advisable to construct the programmed motion 
in the supported phase of a bound, in such a way that it can be per- 
formed with the minimum force or power developed in the log Joints. 

The problems of optimisation of the programmed motion in the supported 
phase and its stabilisation wore investigated in [1] . 

In this work, the problem of performance of the programmed motion 
of the vehicle in the unsupported phase of a bound, constructed in [2] , 
is considered in this work, in the case when various kinds of errors /6 

and perturbations occur. 


The center of mass of the vehicle moves along a ballistic trajectory 
in the flight phase, and its motion is uncontrolled. The angular motion 
of the vehicle housing around the center of mass can be corrected, by 
changing the motion of the legs during flight. For example, by bending 
(unbending) the legs, their moment of inertia changes and, consequently, 
the effect of the transfer of the legs on the angular motion of the 


* Numbers in the margin indicate pagination in the foreign text. 
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housing changes by the law of conservation of angular momentum. 

The possibility of change of the motion of animals and man around the 
center of maos in unsupported motion, by movements of the limbs, has 
quite long been known [3-73 . 

The purpose of the stabilisation of motion of the vehicle in the 
flight phase is to reach the position assigned in the motion program at 
the time of landing, with allowance for the shift caused by deviation 
of the final position of the center of mass of the vehicle from the 
p r q g r amine d position. 

The operation of che navigation and information systems is not 
modeled. It Is considered that, in the supported phase, in the motion 
control system of a bounding vehicle, all current phase coordinates of 
the vehicle are known without error and, in the flight phase, the 
coordinates and angular velocities of the housing, as well as the angles 
and velocities of the leg joints. The control system also knows a model 
of the terrain. At the moment of liftoff from the supporting surface, 
the terrain model in the landing region corresponds to the actual sup- 
porting surface with some degree of accuracy. During flight, the 
terrain model in the landing region is refined, and it becomes completely 
equal to the actual supporting surface immediately before the time of 
landing . 

As to the supporting surfaces, it is assumed that they differ 
little from a section of a plane and that their inclination to the hori- 
zon is small. The terrain as a whole can have a complex form. /7 

Processing of the algorithm for stabilization of the motion of a 
bounding vehicle in the flight phase by mathematical modeling in a 
computer has permitted a quite efficient stabilization algorithm to be 
obtained . 

Formulation of the problem of stabilization of the motion of the 
vehicle In the unsupported phase of a bound and a qualitative analysis 
of various methods of construction of the stabilization algorithm are 
presented in Section 1. The advisability of the use of the principle 
of local determination of supplementary control motion of the legs In 
a forward step in time is shown here. The logic of operation of the 
stabilization algorithm is described in Section 2. The nominal posi- 
tion of the vehicle at the time of landing and the basic characteristics 
of the forthcoming flight phase are determined (its duration, the an- 
gular momentum vector of the vehicle relative to the center of mass, 
etc.). The problem of reorganization of the motion of the vehicle in 
the concluding stage of the flight phase is considered, with account 
taken of more accurate information on the supporting surface obtained 
during flight. At the beginning and end of the flight phase, the 
angles in the leg joints change along the nominal trajectories of 
transfer of the legs during flight, which ensures a shock free liftoff 
and soft placement of the legs on the supporting surface. The angular 
motion of the housing is not stabilized, in this case. During all the 
remaining time, the legs participate at once in two motions: the 

nominal leg transfer motion during flight and the supplementary control 
motions, which provide the change in angular coordinates of the housing 
along the line of transfer, smoothly connecting the corresponding end 


values of the angular coordinates and housing velocity. In Section 3, 
the problem of determination of the supplementary control motion of the 
leg Is reduced to a problem of a special type of quadratic programming, 
the algorithm for solution of which is constructed In Section k. 
Mathematical modeling of the process of stabilisation of the angular 
motion of a bounding vehicle In the flight phase (Section 5) and the /8 
calculation results (Section 6) demonstrate the efficiency of the 
stabilisation algorithm for various conditions of motion of the vehicle 
with perturbations. 

The author thanks D.Ye. Okhotslmskiy for formulation of the 
problem and attention to the work. 

1. Formulation of the Problem 


We will characterize the position of the vehicle housing relative /9 
to the absolute coordinate system (the 0]^ plane Is horizontal), 

by the coordinates of the center of mass of the "'housing £ jri ,£ and the 
angles (Fig. 1): yaw; 0 pitch; y bank. The axes of the Oxyz 

coordinate system connected with the vehicle housing are directed along 
the main axes of inertia of the housing. 



The vehicle has four or six two 
section weighable legs (the number of 
legs Is designated by N). The position 
of the legs relative to the Oxyz axes 
Is determined by the coordinates of the 
points of suspension of the thighs to 
the housing and by the angles (Fig. 2): 
cu is the angle between the Ox axis and 
the plane of the legs; 3^ is the angle 
between the negative direction of the 
Oz axis and the thigh; qi is the angle 
between the thigh and the shank. The 
plane of the legs, formed by the thigh 
and the shank. Is perpendicular to the 
Oxy plane. 

The vehicle moves the Oy axis for- 
ward. 


. . In [2] , an algorithm for the syn- 

’ tas ‘ thesis of the programmed motion of a 

bounding vehicle in the flight phase is considered. From the assigned 
values of the housing coordinates and the absolute coordinates of the 
feet at the times of liftoff from the supporting surface and landing, 
the initial velocity of the center of mass of the housing, the initial 
angular velocity of the housing and the parameters of the leg transfer 
trajectories required to perform the forthcoming flight phase are 
determined. The legs are transferred in the Oxyz relative coordinate 
system, in such a manner that a shock free liftoff and soft placement 
of the legs on the supporting surface are ensured. The initial angular^ 
velocity of the housing compensates the effect of the transfer movement 
of the legs on the movement of the housing around the center of mass 
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during the flight, The programmed motion In 
the flight phase Is complied from the first 
Integrals of the equations of motion (motion /11 
of the center of mass of the vehicle and the 
law of conservation of angular momentum of the 
vehicle relative to the center of mass). 

Actually, In the unsupported phase of a 
bound, the vehicle deviates from the programmed 
motion, as a result of various kinds of per- 
turbation. A variance of the actual dynamic 
and kinematic characteristics of the vehicle 
(dimensions, weight, moments of Inertia) from 
their values In the motion control system Is 
a perturbation . Such a situation arises, for 
Pig. 2. example, If a load Is placed on the vehicle, 

and this is not reported to the motion control 
system. Another type of perturbation involves Inflight errors in 
processing of the leg motions. There are perturbations caused by data 
errors. Up to the time of liftoff from the supporting surface, the 
terrain model in the landing region Is known with a certain degree of 
accuracy. During the flight, Immediately before the time of landing, 
the model of the supporting surface In the landing region Is refined. 

The accurate model of the terrain In the landing region Is used to 
reorganize the motion of the vehicle In the final stage of tho flight 
phase , 

Perturbations in the supported phase by the prior flight phase data 
and time inaccuracy in processing the time of liftoff from the support- 
ing surface result in errors in processing the programmed coordinates 
and velocities of the vehicle at the time of liftoff from the support- 
ing surface. These errors lead to change in the distance of the bound 
and Inflight rolling of the vehicle, and they make it necessary to 
recalculate the characteristics of the forthcoming phase of the flight: 
its duration, angular momentum vector of the vehicle relative to the 
center of mass and the parameters of the leg transfer trajectory, which 
satisfy the conditions of shock free liftoff and the required soft 
placement of the legs on the supporting surface. We will call the in- /12 
flight leg motions, synthesized at the time of liftoff from the sup- 
porting surface, the nominal leg motion. 

The center of mass of the vehicle in the flight phase moves along 
a ballistic trajectory [2]. Therefore, the motion of the center of 
mass is uncontrolled. The motion of the vehjcle around the center of 
mass occurs according to the law of conservation of angular momentum of 
the vehicle relative to the center of mass. In this case, the angular 
motion of the housing can be controlled, by changing the transfer 
motion of the legs during the flight. For example, by bending (unbend- 
ing) the legs, their moment of inertia changes and, consequently, the 
effect of the leg transfer motion on the angular motion of the housing 
changes. We note that, in this case, the angular momentum vector of the 
vehicle remains constant. Only redistribution of the angular momentum 
between the housing and the legs occurs. 
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The purple of stablllautlon of the motion of the bounding vehicle 
in the flight phase Is to reach the position assigned in the programmed 
motion, with allowance for a shift which onuses a deviation of the 
final position of the center of the mass of the vehicle from that pro- 
grammed. The question of determination of the nominal position of th*» 
vehicle at the time of landing will be discussed in detail In flection a. 

Angular motion of the vehicle in the flight phase which satisfies 
these boundary conditions Is provided by Inflight variation of the leg 
motion. 


We will solve the problem of stabilisation of the angular motion 
of a bounding vehicle In the flight phase, by considering the angular 
coordinates and velocities of the housing, as well as the angles and 
velocities In the leg joints to be known without error at the current 
moment of time. 

The motion of the vehicle in the flight phase occurs according to 
the law of conservation of angular momentum of the vehicle relative to 
the center of mass [2] /13 


where & Is the angular velocity of the housing projected on the axes of 
the Oxys relative coordinate system, P=(cii,0i,qq. . .qj*j) are the angles 
of the leg joints, 1 is the matrix for transfer from the ab- 

solute coordinate system to the Oxyz coordinate system connected with 
the housing, TJ' a is the angular momentum vector of the vehicle relative 
to the center of mass In the absolute coordinate system, I and 

Ip are matrices, the elements of which depend only on the angles of 
the leg joints p\ 

We designate the nominal transfer motion of the legs, synthesised 
at the time of liftoff from the supporting surface, p H (t) . We designate 
the angular coordinate vector of the vehicle housing ^=( 140 , 7 ). Let 
t co and be the moments of the start and end of the controlled stage 
of the flight phase, in which t 0 <t c0 <t cl ^t 1 , where t 0 is the moment of 

liftoff from the supporting surface and t^ is the time of landing. 

That motion of the vehicle housing around the center of mass must 
be provided, which satisfies the boundary conditions 


, 9 { % )-% 

YUc,)-% , 


( 1 . 2 ) 


( 1 . 3 ) 


The angles of the leg joints must satisfy the boundary conditions 

^4 / " & (Ye, J f Pi'Ycc . ’ -"4 {£ej, ( 1 . H) 

P ( Ye t J * ,'i lie, J , Ptftt j -k i %}, (1.5) 
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and the phase limitations : 


/V\ 


• *3 ... 

(1.6) 

- H < P 


.■ ft* jj” a * 

(1.7) 

hi itt *.,* • / ** 


’ w ' p ' /; v . 

(1.8) 


Limitations (1.6) -(1.8) involve a given vehicle design. We require 
that condition (1.6) aloe ensure noninterference of the legs of the 
vehicle with each other. 

We note that boundary condition (1.3) on the angular velocity of 

* • 

the housing ^(toi)^^ is satisfied automatically, by virtue of (1.1), 
if boundary conditions (1.2), ( 1 . ‘1 ) , and (1.5) and the first of boundary 
conditions (1.3) are satisfied. 

We will call the difference between the actual and nominal motion 
of the legs the supplementary control motions of the legs 

pft)-fi(t). (i.oi 


We require that the actual leg motion in stabilization of the 
angular motion of the housing differ to the minimum extent from the 
nominal motion of the legs, i.e., 




( 1 . 10 ) 


where u is a certain norm of the vector. For example, 


or 



= m«x /■£ 








£*}>}. 


where ej_ > o,f i >o are weighting factors. /15 

The problem of stabilization of the angular motion of a bounding 
vehicle in the flight phase (l.l)-(l.lO) can be solved as an extreme 
problem with limitations of a general type. Various methods of 
simplification of problem (l.l)-(l.lO) are possible, by means of intro- 
duction of supplementary connections. We consider qualitatively and 
compare some of them. 


Parametric stab I Illation algorithm . 


Vfe assume 


ft a ft ft/ ,4 J, 


0 , 11 } 


l.o., the cuppiomemtax'y control motion of the loco In a fixed function 
of time and parameter vector u (7t~ const) . In this case, fox 1 any 
values of the parameter, there Is 

ft ftft D) a. 'J-'O > “ft <*, 
ft ft ft, $ ) ” & , ft ft ft, , 'ft 1 •’ ft 


l.o., boundary conditions (1,4) -(1,0) ax-e satisfied automatically. 
By substituting (1.11) in (l.b)-(1.8) , wo obtain the region of per- 
missible vaHies of parameters u 


By integrating (1.1) with initial conditions (1.2), over segment 
of time tc[ie 0 ,ft/J , wo obtain the value of the angular coordinates of 

the housing at time t e i, with given values of parameters a 


m,j= /w. 


Then, boundary conditions (1.3) are equivalent to limitations of the 
type, equality to the parameter vector 



(1.13) 


By substituting (1.11) in (1.10), we rewrite the optimising 
functional in the form 


/1 6 




(1.14) 


As a result, we find that the pi'oblem of stabilization of the 
angular motion of a bounding vehicle in the flight phase (1.1) -(1.10) 
is reduced to mathematical programming problem (1.12)-(1.14) . 

For solution of the parametric stabilisation program, a large 
amount of calculation woxlk must be done. This is connected with digital 
integration of differential equations (1.1), with calculation of the 
left sides of limitations of the equality (1.13) type. 

Linear transition stabilisation algorithm . We require that the 
angular coordinates of the housing change along transition lines which 
smoothly connect the boundax^y values of the angular coordinates and 
velocities of the housing (1.2)-(1.3). One method of plotting the transi- 
tion lines is described in Section 2. This approach permits reduction 
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of the dimensionality of extreme problem (,.!)-(!. 10) , which can bo 
solved by determination of the supplementary control motion of the less, 


Local s tabilisation algorithm . We rrqur-'e that the angular coor- 
dinates of the housing also, as In the linear transition stabilisation 
algorithm, change along transition lines. Ta > supplementary control 
motion of the lego Is determined at discrete moments of time T!{(Kn 0 , 1 ,P. . . ) , 
In small forward time steps t J. Wo bourne 


cS - u ({, & ) f m ? 'tli ic L Z, %„], (1.1b) 

where % Is the vector of parameters, which have a constant value In 
a given stop. 

At time T'p ouch values of the parameters must bo found that, with /l' 
taT K+i, angular coordinates and velocities of the housing are on 
the transition linos . 


By substituting (1,15) and the values of the angular coordinates 
and velocities of the housing with the transition lines In angular mo- 
mentum integral (1.1), we obtain a limitation of the type, equality to 
the values of parameters which are written In the following manner 


/ K (l, )~ 0 


(1.1 G) 


Boundary conditions (Ut)-(1.£3) and limitations (1.6)~(1.8) de- 
termine the region of permissible values of the parameters at r <7 

[ h cV K . (1.17) 


Wo require that, at time t k+ ^, the actual motion of the legs 
differ to the minimum extent from the nominal transfer motion of the 
legs (i.e., instead of Integral optimizing functional (1.10), we insert 
a local optimizing functional) 



(1.18) 


By substituting (1.15) In (1.18), we obtain an optimizing functional 
which depends parameters tty 

^ (ef K ) -9tnf, (1.19) 


By solving mathematical programming problem (I.l6)-(l.l8) , we 
determine the values of parameters a K and, consequently, the sup- 
plementary control motion of the legs u(t) with ^ j 


8 



it can be expos tod that the nnallneoo of the atop will permit 
significant niiaplifieation of problem ( 1 .lb)-»i *..18) . 

The basic virtue of tho .local stabilisation algorithm is the lack /18 
of need for digital integration of tho differential equations and, 
consequently, a ve 7 much smaller amount of calculations than In other 
methods of solution of tne problem of stabilisation of the angular 
motion of the vehicle in tho unsupported phase of the bound. 

The most rigid restrictions in selection of tne method of con- 
struction of the algorithm for stabilisation of the angular motion of 
a bounding vehicle In the flight phase is connected with the require- 
ments for high speed of the vehicle motion control computer. The dura- 
tion of the unsupported phase of a bound to a distance for 10 m on the 
surface of the earth is on the order of one and a half seconds. There- 
fore, tho actual requirements for high speed of the computer are made 
only by a local stabilisation type algorithm, 

Qualitative analysis of various methods of construction of the 
algorithm for stabilisation of the angular motion of a bounding vehicle 
In tho flight phase permit it to be concluded that it Is advisable to 
construct a local stabilisation type algorithm. 

2 . Stabilisation bogle 

We designate tho moment of liftoff from the supporting surface 719 
tq and tho moment of landing t^. Tho motion of the vehicle center of 
mass in tho flight phase is uncontrolled. The vehicle center of mass 
moves along a parabolic trajectory [2] 


f* 0 I -*'J * ‘ t* a* 0 '■ 

A n */»<* ,r 


0 


t 


( 2 . 1 ) 


wher§ 3 Ra”(£a> 9 a> Ca) arQ coordinates of the vehicle center of mass, 

O — . 0 

H 7 ,R a are the coordinates and velocity of the vehicle center of mass 
at the time of liftoff from the supporting surface, g is the acceleration 
of gravity. 


Tho purpose of stabilisation is to reach the assigned programmed 
position at the time of landing, with allowance for the shift caused by 
deviation of the final position of the vehicle center of mass from that 
programmed. 


Wo describe the procedure for determination of the nominal position 
of the vehicle at the time of landing. At the moment of liftoff from 
the supporting surface, the model of the terrain In the landing region 
is known with a certain degree of accuracy by the motion control system 
of the bounding vehicle. Let £ s ( n) be the model of the supporting 
surface in the landing region. We require that, at the time of landing, 
the vortical distance from tho supporting surface to the vehicle center 
of mass be the same as in the programmed motion, equal to 


( 2 . 2 ) 
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« ,1 I 1 ■* 

-hero h a °0', A a ,n are ',ho coordinates of the vehicle center of 

noon at the time of landing. 

It follows from the motion of the vehicle center of mans In the 
flight phase that 



vt _ 
(a “ 





Hi - 




4 

*(* 


(a. 3) 


where f n f|_-t 0 Is the duration of the fllgnt phase . 


/SC 


V/e require that the angular coordinates of the housing at the 
time of landing be the same as In the programmed motion. Let 
111 

fr B (g ,n *5 ) be the position of the housing center of mass at the time 
of landing. V/e determine the leg position at the time of landing, by 
using a posture shaping algorithm [13. The supporting contour projected 
on the horizontal plane Q,£p Is a rectangle of fixed dimensions, the 
length of which equals the distance between the suspension points of the 
front and roar lego. The angle between the longitudinal axis of the 
projection of the supporting contour on the horizontal plane and the 
O^p axis (supporting contour orientation) at the time of landing equals 

the angle of yaw of the vehicle housing^ 1 The projection of the 

center of the supporting contour on the horizontal plane lies on a 
line which connects the center of the Initial supporting contour and 


point (s 1 ^ 1 ), at distance ( ( 5 1 ,p 1 ) -«; m i n )/max{l, | | } from point 

where £ nin Is the minimum permissible distance from the 
supporting surface to the housing center of mass, A* Is the tangent of 
the slope angle of the velocity of the vehicle center of mass to the 
horizon at the time of landing. Thus, by using the terrain model, we 
obtain the absolute coordinates of the feet and, consequently, the 
angles of the leg joints at the time of landing, as a function of the 
position of the housing center of mass at the time of landing IT 1 . We 
designate the vector of the direction from the vehicle center of mass 
to the housing center of mass p\ Then, at the time of landing, p -1 is 
a function of K 1 , and the following occurs 


By substituting (2.2)-(2.3) in (2.4), we obtain a system of equations /2I 

relative to PT 1 , in the computer solution of which by the modified 
quickest descent method, we determine the nominal position of the 
vehicle at the time of landing and the duration of the flight phase. 

Because of the conditions of shock free liftoff and softness of 
setting the legs on the supporting surface, the angles of the leg joints 
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change rapidly at the beginning and end of the unsupported phase of 
motion. 'Therefore, in those segments of time, we will not stabilise 
the angular motion of the housing. We designate the times of the start 
and end of operation of the stabilisation algorithm t Q0 and t Ql , with 

V^o^el^r 


At the time of liftoff from the supporting surface, t must de- 
termine the characteristics of the forthcoming flight phase, from the 
phase coordinates of the vehicle at time t , known from the navigation 
system readings * 

The ballistic trajectory of the motion of the vehicle center of 
mass (2.1) Is plotted. The coordinates and velocity of the vehicle 
center of mass at the time of liftoff from the supporting surface are 
calculated easily from the known phase coordinates of the vehicle, 
and they completely determine the motion of the vehicle center of mass 
In the flight phase. 

The nominal position of the vehicle at the time of landing and 
the duration of the flight phase are determined by solution of (2.2)- 
(2.4) . 


The vehicle angular momentum vector relative to the center of 
mass at the time of liftoff from the supporting surface Is calculated. 

The leg transfer trajectory which, ensures shock free liftoff and 
the required softness of setting the legs on the supporting surface with 
low absolute velocity (possibly zero), components of the fixed portion 
of the vehicle center of mass velocity at the time of landing, is 
plotted similar to the way It was done in synthesis of the programmed /22 
motion of the vehicle for the unsupported phase of the bound [2]. 

We will call the Inflight leg motion, obtained at the time of liftoff 
from the supporting surface, the nominal leg motion. 

With ic-[k 0l t e J and r<r/4.4; , the angles of the leg joints 

change along the nominal trajectories, and the angular motion of the 
housing is not stabilized. 

By digitally integrating the equations of the law of conservation 
of angular momentum of the vehicle relative to the center of mass 
(1.1), from the nominal position at the time of landing, with a negative 
time step to t=t c ^ (two steps of integration are completed by the 
Rungs -Kutt fourth order method), the control system determines the 
nominal values of the angular coordinates and velocity of the vehicle 

housing at time t c ]_, which we designate ^ r c ]_i^ r c i (here and subsequently, 

¥=(^,0,y) are the angular coordinates of the housing). 

Let, at time t c0 , there be the following yalues of the angular 

coordinates and velocities of the housing ^co^oo* With t , 

the control system varies the nominal leg motion, in such a way that 
a change in angular coordinates of the housing along the transition 
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Liner, which smoothly connect the boundary value" of 
dl nates and veloc Lt los of the housing (y' ,y., , / and 
ensured. The transition linen are functions of the 


the angular 


type 


ooor- 


then. 


?.* 1 w o 1 »/ ?r* -v; 
« ’ I? {iff. -%) - (t. <-j% > v} 


(a. 10 


where T,«t-t,„. 

L* t_? « 

If, during actual motion at time t=t*, the angular coordinates 
or velocity of the housLng go outside the e corridor of the transition 
Linen 


IY"YJ<£ , (Y~rj*' s > (a. 6) 


a supplementary correction In made, t «t#, and the transition linen 
are recalculated. w 

The Ideas connected with the use of the transition linos are 
consistent with those reported in [8-11] . 

During the flight, before the time of landing, the data system 
refines the model of the terrain in the landing region. By using the 
refined terrain model, at moment t ul , the control system refines the 
nominal position of the vehicle at the time of landing and the dura- 
tion of the flight phase, by solving system (L\2)-(2Jl). The nominal 
leg position Is corrected In the corresponding manner, at IcftuMj. 


In this section, the logic of operation of the algorithm for 
stabilization of the motion of a bounding vehicle In the unsupported 
phase of a bound has been constructed. The problem of determination of 
the supplementary control motions of the legs, which ensure a change in 
the angular coordinates of the housing along the transition lines at 

tc[te 0) h,J , is considered in the next section. 

,j. Problem of Determination of Supplementary Control Motion of Legs 

Let * •%) bo the angles of the leg joints. We 

designate the nominal inflight leg motion, obtained at the time of 
liftoff from the supporting surface, which satisfies the conditions of 
ohocL free lift uff and soft setting, of the legs on the supporting 

surface, by Pj-j ( h ) , and the supplementary control motion of the legs 
by u(t) . We then obtain 


p„ ft) 


rfr ft) cs.D 

We wll.L determine the supp 1 omen t ary control motion of the logs at 
discrete moment s of time t k (K c i?, l ,2. . . ) Ln the forward fLmo stop. Wo 
require that the aeco Loral loan of fho lop, Joint a to constant In the 
forthcoming of op 7e / c\, u< t ] • fa thin case, wo have 


r< (% t , j -• ftz ) t tier* ;^ t v if ft)* 


(3.2) 
( 3 . 3) 


whoro A l K“' i k+ g-'i’,. Lo a otop Ln solution of tho problem of dotermlna- 
flon of the supplementary control motion of tho legs. Tho value of 
Atk Lo selected Ln tho following manner 


A As " nwx /(?•{£*)/ ' 4 

t'°SjJLV 


where e >0 In a small positive value, At raln ,AT max are the assigned 
L Iml to of porml.unl.ble valuer, of Ax^. 

m ♦* 

Ln tho event u(x K ) differs l l.tt.le from u(xK„p) , we find that 
u(i k )At k In a value on the order of or, disregarding which In (3 . a) , 
we obtain 




i f :*f ,■ 


'}* 


>'V. 4 : ./ Af ; 


(3.2*) 


If the anuumption made of the 11 nmoothneuu" of u(t) In not justified, /2b 
the error result Lng from the change to approximate formula (3.2*) will 
Lead to deviation of the angular motion of the vehicle housing from 
the transition linen, and it will become necessary to make a supplementary 
correction. By selection of parameters e,AT jn ^ n and Atj , it can be 

certain that such cases will occur rarely. 

By substituting (3.1), (3.2*), (3*3) , the known nominal transfer 
motion of the legs and the angular coordinates and velocity of tho 
housing, with transition .lines (2.!5) at time x K+1 , in angular momentum 

integral (1.1), wo obtain 


whore 


/i it f‘1)-- i t 
(Pm), 

fPf> I (Pm) - 1 </&)(4 A>. ’ 

Pm ~ A n KM ) m s . 


( 3 . 'I ) 
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Instead of boundary conditions (1.4) -(1.5) and phase limitations (1.6), 
we require satisfaction of more rigid restrictions (3.5) > which auto- 
matically ensures satisfaction of conditions (1.4)-(1,(S) , 


yu/ij FT' 

where U min , U max are fixed quantities, 


(3.5) 


In this case, the coordinate Pj ( j=l,£. . . 3U) transfer trajectory lies 
within the cross hatched region in Pig. 


J * 


We replace restrictions (1.7)- 
(1.8) by the restrictions 


/2 6 



Pig. 3. 


V 

V 


* “ * K 


yn<W 


K 


Him M rr 




i U a!/„ , 


(3.6) 

(3-7) 


where U mln ,U max ,U mln ,U max are 
’ w v s V * w * w 

fixed quantities which ensure satisfac- 
tion of ( 1 . 7) -( 1 • 8) . 


In determination of the supplementary 
control motion of the legs at time t^, 

we will satisfy restrictions (3*5)— (3-7) 
at time It may turn out that. 


at time r K+ ^, the value of u. ( j = l, 2, . . . 3N) proves to be around the 

boundaries of the region of permissible values of Uj (3-5), with a 

velocity directed towards the boundaries and, in the subsequent step, 
limitations (3. 5) -(3 *7) will prove to be incompatible. In order to 
exclude the possibility of the development of such a situation, instead 
of velocity restriction (3.6^, we introduce more rigid restrictions 


-fr — T /-Mat 

V(u) £ M , 


(3.8) 


where 


r#Mtf 


fK"_ , if 
/ / V~*+ VXJA, 


Y I - 


otherwise 


► Wlh 


If u-ju 


fr •**>* / ' , . , 

V (uX"-x 


e v >0 is a fixed positive value. 


otherwise 
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Then , In approaching the boundary of the region of permissible values, 
at distance d u j less than e y , the maximum permissible velocity toward 


the boundarle 
tion to d u j , 


r* 


t 


of the region of permissible values decrease? 
o the velocity along this boundary at d u j=0 . 


In proper 


By substituting the values of u,u (3 • 2)-(3 • 3) In (3 .5)-(3 .8) and 
consolidating these limitations, we obtain 


n 


A m!h 


-1 


ufrji A 




(3.9) 


where 


/f >*** . f r \ 7? — 1 < 7 <=? . 

A “ l [ ~ u kh1^, 

rr~/-htax — ... 7 f , fr >nax J 

iVtiit*)- «(?*>! 3k, % /, 
A"" „„ l [?<?„) if "*•- l fa; 

(V/SL,/ -uKilfoK* 1 '}, 
u’ = D(t.) t Btrj „ 4. 


We require that, at time t K+ ^, the actual motion of the legs 
differ the least from the nominal transfer motion of the legs 

A fa-fa (&»)]** AC #/&)]* J (3.10) 

i-t 

where e i >0,f i >0 are weighting factors. 

By substituting (3.2)~(3.3) in (3.10), we obtain 

ZL { Cc [eft ft K )f * </■ (j/&} (3.11) 

where /2 

p / <?,• 
at - e £ 

As a result, the problem of determination of the supplementary 
control motion of the legs is reduced to quadratic programming problem 
(3.4), (3-9), (3.11), which, by designating x = tr(T K ), we rewrite 
in the form def 
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(3.12) 


3f/ 

ZL Cc { xf X{) , C c >r, 

,e/ 

iN 

£ Q ii' x ; - *! 

$n J ‘ J > 


4 mtM A WAX 

; ■< X.- < A , 


In the event the restrictions of the equality and inequality 
types in problem (3.12) are incompatible, we will look for a solution 
which satisfies the Inequality type restrictions and, with the maximum 
degree of accuracy satisfies the equality type restrictions. For 
this, we solve the problem 


3N n 

*• 

a;u *at : 

S# J 

ZL +2L wy. ->*•#/ 


SH 

L 

i=J 


(3.13) 


where v j>> c ;j_> 0 are weighting factors, y^ are new variables. 

We note that problem (3*13) is practically the same as problem 
(3.12). They are partial cases of the quadratic programming problem 
of a more general type . 


t -s 

JL 4 J{ 2 ; = l 

i-i J J r 

7 mitt _ , — htaX 

4 * 4 ^ 4 

n), I?), 


( 3 . 1*0 


1T1SX 

where Z. are independent variables, Z. =» corresponds to the case 

min 

when Z^ is unbounded above, Z^ =-°° corresponds to the case when Z 
is unbounded below. 


We call the point which satisfies the equality and inequality type 
restrictions the initial approximation of the quadratic programming 
problem. In problem (3.13), the initial approximation is determined 
trivially 

A ‘ 


jr. 


— ‘ v.°= Z 

t Vy lr l J J ' 


ALL the independent variables in problem (i . Id) are bounded above and 
below by finite vaLues. Therefore, in construction of the algorithm 
for determination of the initial approximation of problem (3.1*0* we 

will consider only the ease when all A , max and Z . m * n have a finite 

The surfaces of the optimising functional level of problem (3.1*1), 
are n dimensional ellipsoids, the principal directions of which coincide 
with the directions of change of the independent variables. At the 
same time, the inequality type limitations are limitations on the 
limits of change of the independent variables. These conditions, which 
determine the special form of problem (3.1*0* permit construction of a 
simple algorithm for its solution. 

Thus, the problem of determination of the supplementary control 
motions of the logo, which ensure change of the angular coordinates of 
the housing along the transition lines, is reduced to a special type of 
quadratic programming problem, the algorithm for solution of which is 
constructed in the next section. 

*1. Algorithm for Solution of Quadratic Programming Problem 

An algorithm is constructed In this section, for solution of a 
quadratic programming problem which is a generalisation of problem 
(3.11), (3.13), which develops In determination of the supplementary 
control motion of the legs In the algorithm for stabilisation of the 
angular motion of a bounding vehicle In the flight phase 


* 

c 

at 

h 


£ 



A" £ ■ 


V. * 

it ‘‘ 




(4.1) 


i ( z ft j J ; i fi J 


Algorithms for solution of quadratic programming problems and the 
bibliographies of studies on this question can be found in [12-133. 

The main bulk of calculations In algorithms for solution of quadratic 
programming problems of a general type involve determination of the 
direction of the quickest descent which does not violate the restric- 
tions imposed. The solutions may go beyond the boundaries of the 
region defined by the inequality type of restrictions but, in subsequent 
iteration steps, it may be found that this boundary must be abandoned. 

Problem ( *1 . 1 ) is a special type of quadratic programming problem. 
The surfaces of the optimising functional level 

h 

4—* ~ con&i 

V 
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/33 


are n dimensional ellipsoids, the principal direct lone of which 
coincide with the directions of change of the independent variables* 

At the came time, the inequality type restrictions are reotrletlonc on 
the limits of change of the independent variables . These conditions, 
will ch define the special form of problem ( h . 1) , permit the construct ion 
of a Dimple algorithm for it .a solution. 

Algorithm 1 , solutions of quadratic programming problem (*1.1). 

—P 

We designate Z the solution of problem (4.1), in the case when 

_ p 

there arc no inequality type restrictions. The value of Z is deter- 
mined by the Lagrange indeterminate multiplier method 


% " <?/,[• ( 'll ' "<?<■ ) 

J B * 


(* 1 . 2 ) 


whore X, are Lagrange indeterminate multipliers, which are the solution 
of the J system of linear equations 


L(t A 

{H r°i ^ h ( ; J I't d h; 


(*1.3) 


System (4.3) is nonsingular, if the equality type restrictions are 


linearly independent. 
— P 


rP 


If Z satisfies the inequality type restrictions, 
Z* is the solution of problem (4.1). Otherwise, we take a provisional 

— O 

initial approximation Z , which satisfies the equality and inequality 


type restrictions (the question of the method of determination of Z 

o p 

will be discussed below in this section,) . We connect Z and Z by a 

1 


straight line, and we find point Z between them, in which the 
reaches an inequality type phase restriction in motion from 

_o _P 

Z to Z 


system 


£ 1 =l°y (i p -Z °) , 


(4.4) 


where 






A j 



hint 


, with .?/> f,- 
p <2 


with 
with 


nt<r gP 
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By virtue of (4.4), point Z satisfies the equality type restrictions, 

_P 

since Z and Z satisfy them. 


18 


Wo fix all the variables which reach the Inequality type restric- 
tion.". at point ;T ^ . ’Wo transfer the terms corresponding to them to tin 
equality type restrictions on the right side, recalculating the free 
terms, We exclude the fixed var.lahl.ee from consideration, and we will 
call the remaining variables free. We exclude the linearly dependent 
varlah Leo from the equality type reel riot lone, If they were formed at 
a ronulf of exclusion of the fixed variables. If the number of free 
variables equate the number of equations In the equality t.ype restric- 
tions. 


IT' la the co tut Ion of problem (*'l.l). Otherwise, we cet 0 
.and we repeat everything attain for the free variables. 


tvl 


* 


Theorem l . If there exists Initial approximation , which 
satisfies the equality and inequality type reel r 1 cl lone of problem (*1.1), 
algorithm L, In a fin He, emaller than n number of etepe converter* 
toward the eolation of problem 


Proof . fhe number of variables In problem C h . 1 > equale n. In 
each etep of algor It hm 1, the number of free varlablee deereueee by at 
leaet one. Consequent ly, algorithm 1 etope In a finite, emaller than 
n number of etepe, 


We doe. lunate L the linear multiformity determined by the equality 
type reetr Let tone and si, the region determined by the Inequality type 
reetrtet lone of problem ( h . 1 ) » 


We designate the point towarde which algorithm 1 converges by TT“ 
Wi thout reetr Let Ion of generality, it can be ooneldored that, in tbie 


eaec 


the first n x varlablee reaebed the boundary 


max , 


, l'44V 4 4.4. 4’ K *4^ 4444 4, 444' 4. V 4’ 4 VHV4JV-V. 1444' VWUlHm4,y ( d ^ 83 1 , d , ■ ■ 

and the remaining varlablee lie within region id . Ae wae noted In the 

description of algorithm 1, '£*cL(]Sl, Consequently , In order to 

pi'ove that U# le the solution of problem (*1.1), it is eufftelent to 

show that there do not ex let d irect Lone from point IT#, which do not 
violate the equality and inequality type reetr let lone, along which the 
optimising functional deereaeee . 


.n 


o 


We designate P, (a) ae the projection of vector a in linear multi- 
formity I,. J 

Let variable 2.; t c 0 t n, J reach the boundary of region in the 
£-th etep of the algorithm. Then, 



/• c 


J 


where $ 
At 


le the optimising functional, 
the same time, 

1) *P I . , / ... , f' tinny 


( 4.10 


if 


consequently, 


l i l4L‘J 


( Ip ... MoJ, 


(H. 6) 


For variable:' which do not roach the boundary of region , there in 



( t\a )%+/,,„ n ) . 


(II, 7) 


Lot y bo an arbitrary direction along L, directed from Z A toward 
the interior of region 0. In this 7 ease, wo have 


since the first n Q variables reached the boundary 7 Jio ~Z io max . 
derivative of the optimising functional along y 


The 


(ll. 8 ) 


It follows from (4 *5} =(4.8) that 



(4.9) 


or <j> does not deci’ease along y. 

Since y is an arbitrary direction from Z A , which does not violate 
the imposed restrictions, Z A is the solution of problem (4.1). 

Note . Theorem 1 is not valid for a quadratic programming problem 
of more general form than problem (4.1). For example, we consider the 
quadratic programming problem for two independent variables (Fig. 4). 
There are no equality type restrictions, inequality type restrictions 
are imposed on the independent variables, and they determine the rec- 
tangular region of their change. The surfaces of the optimizing func- 
tional level are ellipses, the principal directions of which do not 
coincide with the directions of change of the independent variables. 

— 0 

In Fig. 4, the initial approximation is designated Z and the point 
towards which algorithm 1 converges, Z - # . The solution of this problem 
is point Z'f’, which does not coincide with Z * . 

We proceed to the algorithm for determination of the initial ap- 

O 

proxlmatlon Z , which satisfies the equality and inequality type re- 
strictions of problem (4.1). We consider the case when the limits 
of change of all the independent variables in the inequality type 
restrictions of problem (4.1) have a finite value. In this ease, 
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without reDtrletlon of generality, It can bo 

considered that z mln nl (1 r 1,2, . . .n) . 

For thia, it Ac sufficient to eubotituto the 
variables 


#\* / a* ***ajr % hurt * *-» 

g. • /? - 4^4£« J 

4 ’ * *1 / * **»»H 

/**' m &\ * 


Further, we will consider that z . max ni, 
a min - 

The inequality type restrictions determine 
region q, which Id a n dimensional cube 


'• *" Ri £ i (t 


(4,10) 


and the equality type restrictions determine the linear multiformity L 


~ (j c 4«? , , kj. 


(4.11) 


Let Z# be the solution of the problem 


% & L, 

il Z* 


(4.12) 


The solution of problem (4.12) Is determined by the Lagrange Indeterminate 
multiplier method, 

Not e . If £*t~ L, ( Z* is the Initial approximation. 


Let 


We prove a series of statements, which is necessary for construc- 
tion of the Initial approximation search algorithm. 

Lemma 1 . Vector IT# is orthogonal to linear multiformity L. 

Proof . It is evident, since Z* Is a tangent of an n dimensional /3b 
sphere with the center at the coordinate origin and linear 

multiformity L. 

By H, we designate the (n-k+1) dimensional linear subspace passing 
through linear multiformity L and the coordinate origin, which is the 



center of n dimensional cube ft . On the strength of Lemma 1, wo have 

fh //.* i**d L t g L <? 1 1 

We designate fi™ the intersection of n dimensional cube ^ anti linear 
subspaae II 

£L h *£L{\H 

Lot; 2 <s H . Then, 5TnX2’H+2‘ j, whore ? % &L . V/e will call |A|the 

distance from point Z to linear multiformity L In subspaae H. flubcpace 
H in divided into two regions by the linear multiformity, The sign of 

A determines to which of these two regions % belongs. By the sign of 

X *» 

A, wo will call those regions II and H . For example, the coordinate 
origin is at a distance of unity from L in H, and it belongs to H”, sine*-' 

We assign ^«MZ*(0<y<l) the point of exit at the boundary of region 
ft, in motion from the coordinate origin to Z* . In this case, one 
variable Z. roaches the boundary of region p , To be definite, we 

JL w 

assume Z S4 io «l. We designate the linear multiformity Z io n l by G 1q . 

Then, is closest to linear multiformity L in H, the bound of 

the section of n dimensional cube ft with linear subspace H. 

Theorem 2 . I £ L 0£L,f4 t L 0 C, ia (\S>^f ¥ (<J> is the empty set). 

Proof . We prove the contrary statement equivalent to this one: /to 

If VI fiSl ~ & . In this case, the bound . V. lies 

entirely on one side of L in subspace H, namely, //", 

Since C.fS--.. is closest to L in H, the bound of the section of cube 

ft with linear subspace II, ft^ also lies on one side of L in H, 

Consequently, LOXb^, By definition, - and -’V/ 

Then, L(l£l >) ,~U}H'\Q.?- Lf]£L . As a result, we obtain Lfl£L"tf. 

Algorithm 2 , determination of the initial approximation of 
problem ( t .1) . 
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Annum© all varialiloo to bo free. 


2. Determine point P, which in the no. lull on of problem C'l .12) 

|| 

for tho free variables. If a '■«!=* , gf is the oolutlon (the stop); 
otherwise, execute paragraph 3* 

3. Find tho point ©f exit on the boundary of n dimensional 

cube 0, for motion from tho coordinate origin to point P, Fix the 
variables reaching the boundary of region ft. Transfer the terms with 
the fixed variables to tho equality type rootx^lctlons on tho right 
side, recalculating the free texrnic. If tho equality typo restrictions 
for tho free variables are Incompatible, tho problem does not have a 
solution (stop). Exclude the linearly dependent equations from the 
equality typo restrictions for the free variables. Pi’oceed to para- 
graph 2 . 


It from theorem 2 that, if the intersection of linear 

multlfoi’mity L and n dimensional cube ft is not empty L 

algorithm 2 converges towards point and, if LflSh^f . 

the algorithm reports the incompatibility of the equality and inequality 
type restrictions of problem (JJ.l). 

5. Stabilisation Process Modeling 

A block diagram for modeling the process of stabilisation of /Hi 

motion of a bounding vehicle in the unsupported phase of the bound Is 
presented in Fig. 5. The modeling was carried out in a BESM-C computer 
in FORTRAN. 

In block 1, from the phase coordinates of the vehicle at the time 
of liftoff from the supporting surface, known from the navigation system 
readings, and the approximate model of the terrain in the landing 
region, obtained from the data system, the characteristics of the 
forthcoming flight phase are calculated. The nominal position of the 
vehicle at the time of landing and the duration of the flight phase 
are determined. Tho ballistic trajectory of uhe motion of the vehicle 
center of mass is plotted. The value of the vehicle angular momentum 
vector and tho angular velocity of the housing at the time of landing 
are determined. The nominal position of the legs in the flight phase, 
which ensure shock free liftoff and soft setting of the legs on the 
supporting surface, is plotted. The times of start t QQ and end t ^ of 

stabilisation of the angular motion of the vehicle housing t 0 < ' b c0 <t c x <1: i 

are calculated, where t ,t^ are the times of liftoff from the supporting 

surface and landing. By Integrating the system of differential equations 
of the law of conservation of angular momentum of the vehicle relative to 
the center of mass, from the nominal position of the vehicle at the 
time of landing, with a negative time step, we find the boundary values 
of the angular coordinates and velocity of the housing at time t 
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Tho otabiliaation algorithm muab cncuro that this pooltion In 
ranched by cho angular coordinators of the heunlng at time t Block 2, 

which Id the algorithm monitor, ceordinatec the operation of ito indi- 
vidual partn , With terto,£< 0 J and U[t< h IJ m the log Joints, the 


nominal inflight motion of the logo Id worked out* With , 

it io required that the angular coordinates of the vehicle housing change 
along the transition lines which smoothly connect the boundary values 
of the angular coordinates and velocity of the housing at times t 

and t Ql - The transition lines are constructed in block 5. She monitor 

chocks the accuracy of motion along the transition lines. If the 
angular coordinates of the housing go outside the e corridor of the 
transition lines, a supplementary correction is prescribed, and the 
transition lines are recalculated. 
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4. 
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landing 


tion of quadratic pro- 

8. 

Input 
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The inflight motion of the legs, which ensures change of the 
angular coordinates of the housing along the transition lines, is 
synthesized in block 4, by addition of the supplementary control motion 
to the nominal leg transfer motion. The control motion of the legs is 


determined locally In step , ahead In time. The problem of 

determination of the supplementary control motion of the legs Is reduced 
to a special type of quadratic programming problem (see Section 3), an 
efficient algorithm for the solution of which Is constructed in Section 4. 

Wo note that the interval between the corrections and the times of 
Intersection of the supplementary control motion of the legs Is at 
least a fixed minimum time segment. 

In block 7, by using more accurate Information on the supporting 
surface obtained by the data system during the flight towards time t cl , 

the control system refines the nominal position of the vehicle at the 
time of landing and the duration of the flight phase. The nominal leg 

motion Is corrected in the corresponding manner, at 

In block 3 , the equations of the mathematical model of the three 
dimensional motion of a bounding vehicle In the flight phase are Inte- 
grated by the Runge-Kutt method. 

6. Calculation Results 


Calculations were carried out for four legged and six legged 
vehicles with the following dynamic and kinematic characteristics. 
Housing mass 1000 kg, housing moments of inertia (In kg*m 2 ) 1=530* 

I =90, I =605* All the legs are the same, and they consist of two 

y z 

thigh and shank sections. Thigh length, 1.2 m, shank, 1.4 m. The 
center of mass of each section is located in the middle. Thigh 2 
weight, 30 kg, shank, 20 kg, moments of inertia of thigh (in kg*ni ): 
I X 2“3.2, I y2 =0.05, I z 2 =3.2; shank: I x x = 3«6, I yl =0.05* I zl =3.6. The 

coordinates of the suspension points of the thigh to the housing of the 


four legged vehicle have the values: X^=0.5 (1=1,2), X^=-0.5 (1=3, 4 1 ), 

Y i =l .25 (1=1,3), Y i =-1 . 25 (1*2,4), Z i =0 (1=1, 2, 3, 4); for the six legged 

vehicle: X^O.5 (1=1, 2, 3), X^-0.5 (1=4, 5, 6), Y ± =1.25 (1*1,4), 

Y i =0 (1=2,5), Y i =-1.25 (1=3,6), Z^O (1=1,2, 3, 4,5,6) . We recall that 

the vehicle moves the 0 axis forward. 

«y 


/4 


The SI unit of measurement system is used. Time is measured in 
seconds, angles in radians. The moments developed in the leg joints are 
measured in thousands of Newton meters. Motion is on the surface of the 
earth, and the acceleration of gravity is 9.8 m/sec 2 . 


The integration step of the vehicle motion model is 0.005 sec. 


The required accuracy of motion along the transition lines, accord 
ing to the angular coordinates and velocity of the housing in the sta- 
bilization algorithm, is 0.01. The minimum permissible time interval 
between two successive corrections and between times of recalculation 
of the supplementary control motion of the legs is 0.025 sec. 
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The Initial and final programmed positions of the yohlole are 
determined by the vehicle posture formation algorithm at the times of 
liftoff from the supporting surface and landing, the operation of which 
was considered In detail In [1]. The projection of the supporting 
contour on the horizontal plane Is a rectangle, the length of which 
equals the distance between the suspension points of the front and back 
logs, and the width of which equals double the lateral extension of the 
legs from the longitudinal axis of the supporting contour r. . The input 
parameters of the posture formation algorithm, the selection of which 
is done, either by the driver of the bounding vehicle, or by the higher 
levels of the motion control system, are: D is the distance of the 

bound (the distance between the projections of the centers of the sup- 
porting contours on the horizontal plane O^gn); I s the direction 

(azimuth) of motion In the flight phase; 


A! 




contour; 
1 


pi. 


is the orientation (azl- 
is the orientation (azi- 
of the angle 


^pl 

Is the tangent 
the center of 


muth) of the Initial supporting 

muth) of the final supporting contour; 
of inclination of the Initial velocity of the center of mass, as well 
as of the model of the supporting surface to the horizon. At the times 
of liftoff from the supporting surface and landing, the vertical dis- 
tance from the supporting surface to the housing center of mass is a 
fixed value ? nom i the pitch and roll angles equal zero (e°=e -y 0 =y 1= °) j 

the yaw angle equals 

( ^°^pl°. The $ ana n 

times are determined from the condition that, 3n the projection on the 
horizontal plane OiCn* these points lie on a line connecting the centers 


the azimuth of the corresponding supporting contour 
The e and n coordinates at the initial and final 


plane 0^£r) * 

of the supporting contours, at distances (S nom “S ln i n )/ ma x{l, 1 
the center of the corresponding supporting contour. Here, C min 
minimum permissible vertical distance from the supporting surface to 


} from 
is the 


the housing center of mass, and A* 
inclination of the velocity of the 
at the corresponding time. 


is the tangent of 
vehicle center of 


the angle of 
mass to the horizon 


For the vehicles under consideration, we assume ? m i n =0i 5 It was 

shown in [1] that, from the point of view of minimization of the maximum 
power developed in the leg joints during the support phase of motion, Ai 

for the four legged vehicle, C nom sl< 7m, r b =*lm are best and for the 

six legged vehicle, ? nom = l ■ 6 m, r^=l.l m. Subsequently, we will consider 

that C nom and r^ have these values. 


As nominal bounds on the horizontal plane we will consider • distance 
D to be 5 and 10 m. The angle of inclination of the initial vehicle 
center of mass velocity to the horizon is 45° • The angle of rotation 

1 0 

of the vehicle in the flight phase piA pj_ is 0 and 0.3 radians. 

The direction of motion of the vehicle in the flight phase is \p ~= 

° 1 * 
l/2(i|jpi +i/> ^ ) • We select the absolute coordinate system in such a 

way that £°=0, n°=0, 1 4p 1 o=5 0. 
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Wo Investigate the efficiency of processing of various types of 
perturbations for these modes of motion. The required accuracy of 
motion along the transition lines In the algorithm for stabilisation of 
the angular motion of the housing In the flight phase Is 0.01 radian. 

We estimate the efficiency of the stabilization algorithm with pertur- 
bations, from the accuracy of processing of the nominal angular coordi- 

111 

nates of the housing at the time of landing. We designate 6 ip ,60 ,<$y 
as the deviation of the actual angular coordinates of the housing from 
nominal at the time of landing. The stabilization algorithm efficiently 
processes inflight perturbations, If max{ | 1 ( , | ; 5 y 1 ( }^0 .01, i.e., 

the deviation of the actual angular coordinates Oi he housing from 
nominal at the time of landing does not exceed the required accuracy of 
motion of the angular coordinates of the housing along the transition 
lines. We will consider 0.03 radian (approximately 1.7°) accuracy of 
processing the nominal angular coordinates of the housing at the time 
of landing to be satisfactory. 


Let the values of the dynamic characteristics (masses and moments 
of inertia) of the leg sections of the vehicle be loaded into the 
motion control system with errors. We consider the effect of pertur- 
bations due to these errors. Let m^,!^ be the actual masses and 

moments of Inertia of the j -th section of the i-th leg (j=2 corresponds 
to the thigh, j=l, to the shank), and be the masses and 

moments of Inertia of this section loaded into the control system. We 
consider the case when the magnitudes of the relative errors are the 
same with respect to the masses and moments of inertia of all sections 
of the legs and equal to 


/*»7 
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W If V 



V 






For the modes of motion listed above, errors of + 5 % ( [ Sm^J^O . 05) are 

processed efficiently. In this case, the accuracy of processing the 
nominal angular coordinates of the housing at the time of landing, 
practically, is not over 0.01 radian. With relative error fim^ of -15% 

to 25 % , the accuracy of processing the nominal angular coordinates of 
the housing at the time of landing is not over 0.03 radian. The 
maximum error arises in pitch angle 60 1 , and the nominal values of the 
bank and yaw angles are processed practically ideally. Errors <Sm H >0 


are processed better than negative errors. In other words, it is 
better if the motion control system considers that the legs are heavier 
than actual, and not the contrary. These results are valid for the 
four legged and six legged vehicles. The accuracy of processing the 
nominal angular coordinates of the housing at the time of landing vs. 
fim^, for a bound to a distance of 10 m and an angle of rotation in the 


flight phase A\J; 


pi 


=0.3 radian is shown in Table 


1 . 


We consider another type of perturbation, connected with the in- 
accuracy of processing the leg motion in the flight phase. It is 
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TABLE 1 


a EOHf»t!Hftn ""'Hl'Kn 

f • . . • TTk-ir 'Musi aauaxai “ ’ * _ o S=?T norm: Mnanav 

m h , ..... -™- — . — — ~ -=m ! 7TT, 1 7~. 


! 

se‘ h’H* 


E0 l j S/‘ j 

: !’ ’ “ 

i c.05 *; * 

0.075 !‘ s 

> C.I0 j K 
, 0.125 » 

• 0,15 V » 

' 0.175 * 

0.20 ! * 

! C.25 i; x 
C.3C * 

X X 

x i x 

-0.012 ■ X 

-0,013 X 

-0.015 f X 
-0.CI5 X 

-0,019 • X 
-0.C20 X 

-0.032 X- 

x 

X 

: 

s 

X 

X- 

X 

i 

X I * 1 

X X i 

-0.0II i x, 1 

-0.014 X ! 

-0.016 x j 

-0.020 x ; 

-0,021 1 * 

-0.021 ! x 

-0.027 x 

4 

1 — 0 ,05 j X 

! -0.075 * 

• -0.10 , X 

. -C.I25 X 

: -0.15_ ;• x 

i -0.175 X 

i -0.20 • X 

s ! 

X * 

0.015 X 

0,023 . X- 
0.028 ; -0.CXI 
0.025 ! -0.013 
0.033 : -0.013 
j 0.070 ; -0.02'* 

! x 

i K 

1 x 

X 

X 

X 

X 

X 

0.016 

0,023 

0.027 

0.024 

0,027 

0.032 



J 

-c*oic 

-0.0IC 

-0.010 

-0.0II 


* - error less than 0.01 


Key: a. Pinal error 

b. *1 legged vehicle 

c. 6 legged vehicle 
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considered that, at the er , ror i n processing the required /H 9 

accelerations in the leg joints is an unbiased normal distribution of 
a random quantity. We recall that the minimum permissible step of 

recalculation of the leg motion is 0.025 sec. Let P ^ n ( t ) (j=l, 2, . . .3N) , 

where . . . ,q^, be the required acceleration of the leg joints 

at time t. We set the actual accelerations of the leg joints equal to 

p(t)=p n (t ) (1-i-v ) , where v is the normal distribution of a random quantity 
with a mathematical expectation of zero and a root mean deviation of 
a. The calculation results show that, for the modes of motion listed 
above, an unbiased random error of the required acceleration of the 
leg joints, on the order of their nominal value (with a<l) , has prac- 
tically no effect on the accuracy of processing the nominal angular 
coordinates of the housing at the time of landing max{|Sifj ■*•[ , | 6 0 1 [ , | 6y^| } 

40 . 01 . 


At the time of liftoff from the supporting surface, information on 
the supporting surface in the landing region is known with a certain 
degree of accuracy to the motion control system of the bounding vehicle. 
During the flight, the model of the terrain in the landing region is 
refined. The precise model of the terrain is used to reorganize the 
vehicle motion in the concluding stage of the flight phase, at 

We investigate the effect of data errors on the accuracy of processing 
the nominal angular coordinates of the housing at the time of landing. 

We consider bounds on a horizontal plane without turning in the flight 
phase (A^p^=0), to distance D between 5 and 15 meters, to be nominal 

motion. The angle of inclination of the initial vehicle center of mass 
velocity to the horizon is 45° • During the flight, the motion control 
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system determines that the altitude of the supporting area differs from 
the nominal value hy the amount <$h. A data error on the supporting 
area altitude 6h for these modes of motion lead only to error In /52 

processing the nominal pitch angle at the time of landing (SO 1 , and the 
bank and yaw angles are processed practically Ideally. The areas of 
values of data errors fih, in which 1 50 - 1 1 does not exceed 0.01, 0.02 and 
0,03 radian are shown in Pig. 6 , for a four legged vehicle. The per- 
missible errors <5h depend essentially on bound distance D. We designate 
Dp the distance between the front legs of the vehicle at the time of 

liftoff from the supporting surface and the back legs at the time of 
landing. For the vehicles considered In this section, D^=D-2.S m. For 

the four legged vehicle, we obtain a value of | 5 0 ^ { <0 .01, with 
-0 .0l6<<Sh/D f <0 .028, i.e., such data errors have practically no effect 

on the accuracy of processing of the nominal angular coordinates of the 
housing at the time of landing. At -0 .065<6h/D^<0 .08, there, is satis- 
factory accuracy in processing the nominal angular coordinates of the 
housing at the time of landing, |<S6 1 |<0.03* Similar results, obtained 
for a six legged vehicle, are presented in Fig. 7. Here, the permissible 
data errors Sh are somewhat less than for the four legged vehicle. This 
is because the total weight of the legs is 1.5 times greater in the six 
legged vehicle and, consequently, the errors of processing the nominal 
angular coordinates of the housing at the time of landing, due to re- 
organization of leg motion of the vehicle in the concluding stage of 
the flight phase, are greater. 



Fig. 6 

To bring about the supported phase of the motion of a bounding 
vehicle, in this case, the active perturbations cause errors in 
processing the programmed coordinates and velocity of the housing at the 
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Pig. 7. 


time of liftoff from the supporting surface. We Investigate the ques- 
tion of the effect of these errors on the accuracy of processing the 
nominal angular coordinates of the housing at the time of landing. /54 

We consider a bound on the horizontal plane to a distance of 10 m as 
an unperturbed motion. The angle of rotation of the vehicle In the 
flight phase is 0.3 radian, the angle of inclination of the initial 
vehicle center of mass velocity to the horizontal is 45° • We will 
assign an error of one angular coordinate or the velocity of the hous- 
ing at the time of liftoff from the supporting surface and determine 
the accuracy of processing the nominal angular coordinates of the 
housing at the time of landing, by the stabilization algorithm. The 
calculation results for four legged and six legged vehicles are 
presented in Table 2. In the event the initial errors of the angular 
coordinates of the housing do not exceed modulus 0.05 radian, and the 
initial errors of angular velocities of the housing do not exceed 
modulus 0.1 rad/sec angle of pitch or 0.15 rad/sec angles of yaw and 
bank, the stabilization algorithm processes the nominal angular coor- 
dinates and velocity of the housing at the time of landing to within 
0 .01 radian. 


We consider in detail the process of stabilization of the angular 
motion of a bounding vehicle in the unsupported phase of the bound, in 
the case when the unperturbed motion is a bound to a distance of 7-5 m 
in the horizontal plane. The final supporting area is turned by 0.2 
radian relative to the initial one, and the deviation of the direction 
of motion of the vehicle from the azimuth of the initial supporting 
area is 0.1 radian. The bound is made at an angle of 45° to the 
horizon. A 100 kg load (10# of the weight of the housing) is placed 
on the vehicle. The motion control system does not know of its presence. 
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TABLE 2 
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* - error less than 0.01 


Key: a. Pinal error c. 6 legged vehicle 

b. 4 legged vehicle d. Initial error 

In this case, the moments of inertia of the housing change in the 

following manner (in kg'm 2 ): ^I =25, ^1 = 5, <$I =2 5- The perturbations 

x y z 

introduced by the load cause deviations from the transition line and 
make it necessary to make supplementary corrections, both in the sup- 
ported phase of motion, and in the flight phase. The vehicle coor- /59 

dinates and velocities at the time of liftoff from the supporting 
surface were obtained by the operation of the motion control system 
in the supported phase. The perturbations introduced by the load in 
the supported phase caused inaccuracy in processing the programmed 
coordinates and velocities of the housing at the time of liftoff from 
the supporting surface by the algorithm for stabilization of the motion 
of a bounding vehicle in the supported phase. The process of stabiliza- 
tion of the motion of a bounding vehicle in the flight phase is shown 
in Pigs. 8-11. The time dependence of the angular coordinates of the 
vehicle housing are presented in Pig. 8. The stabilization algorithm, 
efficiently processes the perturbations, and the angular coordinates 
of the housing at the time of landing scarcely differ from the nominal 
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values. The angle.’ of the leg joints and their 
velocities change continuously during motion 
of the vehicle (Pig. 9-10) • The time depend- 
ence of the forces developed in the leg joints 
(in thousands of Newton meters) Is presented In 
Pig. 11. At the start and end of the flight 
phase, the angular motion of the housing is 
unstabilized, and the legs participate only In 
the basic transfer motion of the legs (see 
Section 2). The maximum moments developed In 
the leg joints in these time segments are 
scarcely less than the maximum moments developed 
in the leg joints during the flight. Con- 
sequently, the necessity for completing the 
bound in the unsupported phase, with the basic 
transfer motion of the leg, the additional 
control motion of the legs does not result In 
an Increase in the forces developed in the leg 
joints. The moments developed In the leg 
joints in the supported phase of motion were 
obtained in [1] . In the flight phase, the 
forces developed in the leg joints are approx- 
imately 5 times less than in the supported 
phase of motion. 

By using the similarity and scale methods 
in the mechanics [14] , the results obtained can 
be generalized. With a proportional & times /6 0 

change of all linear dimensions and a m times 
change in the weights of all the components of 
the bounding vehicle, the basic characteristics 
of motion of the vehicle change in the follow- 
ing manner: 

flight phase duration, /3T times ; 

__ angular velocity of housing and velocities of leg joints, 

1//3T times ; 

accelerations of leg joints, 1/Z times; 
moments in leg joints, m& times. 

For example, with a twofold decrease in linear dimensions, let the 
mass decrease in proportion to volume (the vehicle has the same 
specific weight). Then, the moments in the leg joints decrease 16 times. 

The calculation results show the efficiency of the algorithm for 
stabilization of the angular motion of a bounding vehicle in the flight 
phase, for various modes of motion with perturbations. 

Conclusion 






Fig. 8. 


The problem of stabilization of the angular motion of a bounding /61 
vehicle in the flight phase is investigated in the work. Control of 
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Fig. 11. 

the angular motion of the housing In the unsupporteo phase of the 
bound Is accomplished by a supplemental control motion of the legs. 

The following conclusions can be drawn from the results of the work: 

1. the advisability of the use of the principle of local 
determination of the supplementary control motion of the legs In the 
stabilization algorithm was demonstrated; 

2. an algorithm for stabilization of the angular motion of 
a bounding vehicle in the flight phase was constructed and utilized by 
computer; 


3. the results of mathematical modeling of the process of 
stabilization of the motion of a bounding vehicle in the flight phase 
show the efficiency of the stabilization algorithm, for various modes 
of motion of the vehicle with perturbations; 

4. the magnitudes of the moments in the leg joints, devel- 
oped in the flight phase, are approximately 5 times less than in the 
supported phase of motion; 

5. the necessity for making inflight supplementary control 
motions of the legs does not lead to an increase in the maximum moments 
in the leg joints. 
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